Observation of Intense Second Harmonic Generation from M0S2 Atomic Crystals 
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The nonlinear optical properties of few-layer M0S2 two-dimensional crystals are studied using 
femtosecond laser pulses. We observed highly efficient second harmonic generation from the odd- 
layer crystals, which shows a polarization intensity dependence that directly reveals the underlying 
symmetry and orientation of the crystal. Additionally, the measured second-order susceptibility 
spectra provide information about the electronic structure of the material. Our results open up new 
opportunities for studying the non-linear optical properties in these novel 2D crystals. 
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The family of 2-dimensional (2D) transition metal 
dichalcogenide semiconductors MX2 (M for Mo or W, 
and X for S, Se or Te) has attracted much recent at- 
tention. They are direct band gap semiconductors that 
possess properties very different from graphene [THTT]. 
For instance, single layer M0S2 shows efficient light emis- 
sion [5, 6 , optical control of valley polarization [THS], en- 
hanced many-body interactions [12] , and transistors with 
high on-off ratios [I3j [14] . These properties are not only 
important for applications in electronics and optoelec- 
tronics with 2D crystals, the associated valley-dependent 
physics and topological transport phenomena [15] are 
also of fundamental interests. 

The lack on an inversion center in the crystalline struc- 
ture of the material (so that energy gaps at the K and 
K valleys with finite but opposite Berry curvatures are 
introduced) is the underlying reason for developing inter- 
esting properties. Actually, the inversion asymmetry of 
the material not just affects its electronic and linear opti- 
cal properties, but also gives rise to a finite second-order 
optical nonlinearity (x (2) + 0) [T6| fT7] . The finite x (2) 
provides opportunities for new applications in optoelec- 
tronics devices (e.g. coherent control of valley-polarized 
currents) that are not possible with graphene. 

In this Letter we report the observation of intense sec- 
ond harmonic generation (SHG) from odd-layered M0S2 
thin films. In contrast, very low SHG is observed for sam- 
ples with even-layer numbers due to the restoration of in- 
version symmetry. Moreover, the 3-fold rotation symme- 
try of the crystal produces a characteristic 6-fold polar- 
ization dependence pattern for the SHG intensity, from 
which we have developed an optical method for imaging 
the underlying crystalline structure and domain orienta- 
tions with sub-micron resolution. We have also measured 
the x^ spectra by varying the incident pump photon 
energy. Comparison to absorption measurements shows 
that the observed efficient SHG comes from the reso- 
nantly enhanced virtual optical transitions in the near 
UV range. 

The M0S2 mono and few-layers were obtained by 
mechanical exfoliation of natural bulk hexagonal M0S2 
which is deposit in Si covered with 300 nm Si02 sub- 



strates [1 or on transparent amorphous quartz (SPI Sup- 
plies Inc.). We have characterized the M0S2 layer num- 
ber with a combination of optical microscopy, AFM [18] , 
Raman spectroscopy [18] and photoluminescence imag- 
ing [5]. For the second harmonic imaging we used a 
140 fs Ti-Sapphire oscillator (Coherent Chameleon) with 
80 MHz repetition and tunable wavelength from 680 nm 
to 1080 nm which is directed to a confocal scanning laser 
microscope (Olympus FV300) modified for two photon 
excitation. The laser is focused on the sample at nor- 
mal incidence by a 60 x objective (NA 0.75). The back 
reflected signal is then directed to a dichroic mirror and 
a thin band pass centered at the SH wavelength to com- 
pletely remove the laser scattered light and the SH signal 
is detect by a photomultiplier tube. 

Fig Ilia) shows the optical images of mono- and few- 
layer M0S2 samples, where the colors indicate the number 
of layers [15] . Fig.[IJb) shows the sample characterization 
by atomic force microscopy (AFM) at the dashed trian- 
gle area marked in Fig. pja). By measuring the height 
of M0S2 relative to the substrate [18 we found regions 
of M0S2 monolayer (1L), bilayer (2L) and trilayer (3L) 
respectively, as indicated by the labels at Fig. []Jb), see 
supplemental materials for details. Fig.JIJc) presents the 
sample second harmonic image, same area as in fig. Ilia), 
showing well defined regions where the second harmonic 
intensity varies. To unambiguously show that the strong 
optical emission in Fig. [IJc) is due to the optical second 
harmonic generation, we measured its intensity depen- 
dence on pump laser power as shown in Fig. []Jd). Since 
SHG is a second order nonlinear optical process, the SHG 
intensity scales quadratically with the fundamental pump 
intensity [I6j [17 . This is in agreement with our results 
shown in Fig. [lid) where the log scale plot shows a linear 
dependence with a slope equals to 2.00 ± 0.02. 

To visualize the spatial variation of the second har- 
monic intensity as a function of the M0S2 ultra-thin film 
position we plot the SHG intensity profile in log scale, 
Fig. file) at the marked yellow line in Fig. [Tic). It is 
intriguing to observe that a stronger signal is present 
for the monolayer M0S2 when compared to the trilayer 
M0S2, for the pump wavelength at 800 nm, considering 
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FIG. 1. (color online) (a) Optical microscope image of the M0S2 thin film, (b) AFM image of the same M0S2 thin film made 
within the dashed triangle shown in part a, scale bar 1 /mi. By measuring the height relative to the substrate, the number of 
M0S2 layers are determined as 1L for monolayer, 2L for bilayer and 3L for trilayer M0S2. (c) Second Harmonic image collected 
from the same M0S2 thin film, pump laser wavelength at 800 nm, 1.55 eV. Brighter colors means stronger SHG intensity, and 
the scale bar is 5 /xm. (d) Log vs. Log plot of the power dependence of SHG from the X marked position at part c. The circles 
are the measured data and the red line is a linear fit of the data with extracted slope of 2.00 dz 0.02. (e) Intensity profile of the 
SHG image from left to right at the yellow line shown in part c. (f) Atomic structure of monolayer (left) and bilayer (right) 
M0S2. 



that less material is present in the former case. Also 
the bilayer M0S2 shows very low second harmonic gen- 
eration, almost two orders of magnitude lower than the 
monolayer. The small observed signal may come from 
the boundary between the layers [20j[2T]. The x^ tensor 
holds the information of the crystallographic properties 
of the crystal [16] [17] and for SHG without phase match- 
ing determines the SHG intensity. This is true for the 
present case where the sample thickness is much smaller 
than the light wavelength. Hence, in order to understand 
this spatial intensity variation, we should determine each 
element of the tensor to directly probe the crystal struc- 
ture of two-dimensional M0S2. 

The bulk M0S2 crystal has a trigonal prismatic struc- 
ture with Bernal stacking [22], thus it belongs to the 
P6s/mmc non-symmorphic space group [22] [23] , with 
an inversion symmetry operation in the middle of the 
two M0S2 monolayers. The monolayer and bilayer M0S2 
structures are shown at Fig. [ljf ) . As the layer number 
decreases down to one monolayer of M0S2, it is conve- 
nient to make the point group representation of the unit 
cell, which is comprised by three atoms (one Mo and two 
S). Therefore we assign the symmetry of the monolayer as 
D 3 h, since the inversion symmetry present in bulk M0S2 
is lost. However, in bilayer M0S2 the inversion symmetry 
operation is again present, nonetheless due to the lack of 
translation symmetry along the z-axis, the point group 
is D3d. In general, even number of layers belongs to D3^ 
point group and odd number of layers to the D3^ point 



group (similar to few-layers graphene point groups [24]). 
This simple approach can explain qualitatively the ab- 
sence of second harmonic generation in bilayer M0S2: due 
to the presence of inversion symmetry the x^ tensor is 
equal to zero, thus no second harmonic is generated for 
bilayer M0S2 within the dipole approximation [I6j [17] . 
Generally speaking, there is no second harmonic genera- 
tion for even M0S2 layer number, but it can be generated 
in odd number of M0S2 layers. We should note that, we 
also studied monolayer graphene exfoliated on transpar- 
ent amorphous quartz and observed no significant SHG 
signal, as is expected from the symmetry of the crystal. 

Due to the sensitivity on crystal symmetry, 
polarization-resolved SHG measurements provide 
important crystallographic information of M0S2 atomic 
layers. To measure the polarization dependence of the 
SHG, we have placed the sample on top of a precision 
rotation stage, which is placed under the microscope 
objective. The laser (at 800 nm) is focused on the sample 
at normal incidence and with a fixed linear polarization 
that is in the plane of the sample. The sample is initially 
with the crystallographic axis at an arbitrary angle to 
the incident polarization. We have measured one SH 
image for each sample angle, at 5 degree steps, for the 
SH polarization parallel to that of the incident laser (a 
polarization analyzer is placed before the detector to 
select the SH polarization). 

Figure [2Va) shows the experimental results of the SHG 
polarization dependence for the monolayer of M0S2 in 



Fig. [TJc), where the second harmonic intensity is plotted 
as a function of the sample rotation angle. A clear 6- fold 
pattern can be observed at Fig. [5|a). In figure [2^b) is 
a sketch of a top view showing the angle of the M0S2 
crystallographic orientation to the incident laser polar- 
ization direction (e w ). Figure, ^[c) shows the SH images 
for the monolayer and trilayer M0S2 samples together 
with the crystallographic orientation of the crystal lat- 
tice obtained from the polarization data in Fig.|2Va). To 
obtain the crystallographic orientation, we first describe 
the electric field of the generated second harmonic light 
(E(2cj)) along a given direction (e2 W ) in terms of the x^ 
tensor and input light polarization vector (e w ) as p!6l [T7]: 



E(2cj) • e 2o j = Ce 2o 
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where uj is the laser frequency, 2u is the SH frequency and 
C is a proportionality constant which contains local field 
factors determined by the local dielectric environment. 
Equation [I] fully describes the polarization dependence 
of the SHG, and its specific form will depend on the x^ 
tensor. For the odd-layer with D3^ point group symme- 
try, the second order susceptibility tensor has a single non 

zero element PHE]: XMoS 2 ^XxL=-Xxyy=-Xy%=-Xyxy, 
where x corresponds to the armchair direction, since it 
has a mirror plane symmetry p~6j [17] , and y is the zigzag 
direction, see Fig. [2}d. Thus, using Eq. [I] we obtain the 
resulting dependence of the generated second harmonic 
electric field as a function of the sample angle for a pump 
laser polarization (e w ) parallel to the analyzer (e^), that 
can be expressed as p~6j [17] : 
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where <p is the angle between input the laser polariza- 
tion and the x direction, and <fio is the initial crystallo- 
graphic orientation of M0S2 sample. It is worth men- 
tioning that for the analyzer at crossed polarization with 
the input laser the second harmonic electric field is pro- 
portional to sin(3(j) + 0o)j hence the total SH intensity 
(without analyzer) is constant with 6. The intensity of 
the generated second harmonic light as a function of sam- 
ple angle is obtained from Eq. [2] [16] H7j . resulting in 
I20J oc cos 2 (30 + 4>q). Using this intensity dependence of 
the second harmonic with sample angle, we have fitted 
the experimental data as shown in Fig. [2|a). The fitting 
shown by the red lines agrees with the experimental data 
within 0o = 16 ± 5 degrees, indicating that the armchair 
(x) of M0S2 was initially rotated with respect to input 
laser polarization. We should note that as our measure- 
ment is phase insensitive, there is an arbitrariness of 60 
degrees in the definition of the x-axis. 

This result shows that from this simple approach, 
one can measure the crystallographic direction of M0S2 
purely by optical means. As an example of the appli- 
cation of this technique, we have shown in Fig. p[c) the 
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FIG. 2. (color online) (a) Polar plot of the second harmonic 
intensity from monolayer M0S2 as a function of the sample 
angle, (b) Top view of the M0S2 crystallographic orientation 
with respect to the incident laser polarization (e w ). (c) SH 
image of the same sample as in Fig. [it showing the crystallo- 
graphic direction of the monolayer and trilayer M0S2 sample 
determined by our polarization measurements. 



SH image of the monolayer and trilayer M0S2 sample to- 
gether with the crystallographic orientation of the crys- 
tal lattice determined by the polarization measurement 
in Fig. |2|a). This method can be used for future appli- 
cations in materials where it is important to determine 
the crystallographic directions. Two well known exam- 
ples are the 2D heterostructures (i.e. stacks of different 
2D crystals) and nanoribbons [25-27 where the physical 
properties depends on the relative 2D crystal stacking 
orientation and the edge atomic structure, respectively. 
The SHG can provide also valuable information about 
the electronic structure of crystals. Hence, we have mea- 
sured SHG intensity as a function of the incident laser 
wavelength for the monolayer and trilayer. For these 
measurements the laser beam was redirected to a home- 
built microscope with 60 x objective, 0.75 NA, giving a 
spot size of 1.2 fivn in diameter at the sample focal plane. 
The back scattered light is collected and sent to a spec- 
trometer and detected by a charge coupled device (CCD) 
with a set of filters to remove the pump laser light. In 
order to avoid any spectral dependence of the SHG signal 



coming from our experimental setup, we have used the 
same laser power and we have verified that the spot size 
does not change significantly for the different laser lines 
used. Also, the spectral dependence of the spectrometer, 
collection optics and CCD camera were normalized by 
using a calibration lamp (Ocean Optics). 

Figure [3ja) shows an optical image of the M0S2 sample 
on transparent amorphous quartz where we have char- 
acterized the mono-, bi- and tri-layers by AFM. Fig- 
ures [3Fb)-(d) are SH images of such sample taken un- 
der different incident photon energies where is possible 
to observe that by decreasing the laser energy the mono- 
layer and the trilayer SH intensity becomes very similar, 
Fig. [3^c), and at Fig. J3jd) the monolayer intensity be- 
comes slightly smaller as compared to the trilayer. Fig- 
ure [3le) shows the SH signal for the monolayer and tri- 
layer samples as a function of the SH energy, where res- 
onance enhancement peaks are clearly observed. To ex- 
tract this spectral dependence we normalized the SHG 
intensity from M0S2 (ImoS 2 ) t° that from alpha-quartz 
{I quartz)- For samples on transparent amorphous quartz 
substrates with bulk inversion symmetry, the susceptibil- 
ity ratio can be written as [20j |21] (details in supplemen- 
tal materials): 



,( 2 ) 



*MoS 2 _ 7T(N.A.) 2 f y/l M oS 2 (2u)/lMoS 2 (w) 
Xquartz " Ank ^ Nd V Quartz (2uj) / 1 quartz (uj) 



(3) 



here n is the refractive index of quartz (note that we 
have neglected the slight spectral variation of n and 
the small difference of n between amorphous and crys- 
talline quartz), / is a numerical constant that depends 
on the numerical aperture of the objective (N.A.), N is 
the number of layers, d is the M0S2 interplanar distance 

(0.67 nm), ko=27r/\ is the magnitude of the vacuum fun- 

(2) 
damental wavevector and Xquartz * s the bulk susceptibil- 
ity of alpha-quartz. Fot the bulk quartz, we have used 
the expression for back reflected SHG intensity from the 
alpha-quartz surface [28] . 

The X { MoS 2 /Xq 2 u\rtz for both the monolayer and tri- 
layer samples are plotted in figure [3le) as a function of 
the SH energy, frequency 2uj. Also shown is the normal- 
ized linear absorption spectra of the corresponding sam- 
ples (details of measurement setup in the Supplemental 
Materials). We can see that over the specified spectral 
range, there are three peaks known as the A, B and C 
excitons. They originate from the direct optical tran- 
sitions at the K (for the A and B features) and the T 
(for the C feature) point of the Brillouin zone [5]. A 
significant enhancement of Xm s 2 * s observed with spec- 
tral positions very well matched with the C-resonances 
in the linear absorption spectra (~ 2.8 eV for the mono- 
layer and ~ 2.7 eV for the trilayer sample). It suggests 

(2) 
that the resonance enhancement of XmoS * s due ^° ^ ne 

increased density of electronic states at the V point, see 
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FIG. 3. (color online) (a) Optical image of the few-layer 
M0S2 sample on transparent quartz substrate. The number 
of layers are indicated as measured by AFM. (b) SH image 
for the same sample with pump laser wavelength at 860 nm, 
1.44 eV, (c) 900 nm, 1.38 eV, and (d) 960 nm, 1.29 eV. (e) 
Second order susceptibility, left scale, of a sheet of M0S2 for 
monolayer (blue circles) and trilayer (green circles) samples 
with respect to the bulk quartz as a function of pump laser 
energy. Measured linear absorption (right scale) spectra of 
monolayer (blue line) and trilayer (green line), with the la- 
beled optical transitions. The inset shows the diagram of 
the second harmonic enhancement in monolayer M0S2 where 
the two-photon energy is in resonance with the C absorption 
peak. 



inset of Fig. |3Fe). The different resonance position for 

the monolayer and the trilayers is thus a consequence of 

the different energy gaps at the T point of the Brillouin 

zone, that is a modification of the electronic structure 

from quantum confinement at reduced sample thickness. 

(2) 
The SHG from the trilayer Xm s 2 * s consistently smaller 

than that of the monolayer. Away from the resonances, 
the XmoS °f the- two odd layers become similar. Such 
an observation demonstrates the perfect cancellation of 
second harmonic (SH) dipole moments in an AB-stacked 
Bernal unit cell. 

In summary, our demonstration of strong SHG from 
monolayer and trilayer M0S2, although being only a few 
atoms thick, shows its potential use for nonlinear opti- 
cal device applications. Also the SHG dependence on 
the crystallographic axis and electronic resonance effects 
open up an optical way to characterize the crystal struc- 
ture of the many novel 2D crystals. Furthermore, the 
ability to distinguish different crystallographic axis can 
be applied to improve materials properties, such as stack- 
ing different 2D materials with atomic registry and edge- 
defined nanoribbons. 
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